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We derive the spectral density of a sinusoidal carrier phase modulated 
by a random baseband pulse train in which the signaling pulse duration 
is finite and the signaling pulses may have different shapes. The spectral 
density is expressed as a compact Hermitian form in which the Hermitian 
matrix is a function of only the symbol probability distribution, and the 
associated column vector is a function of only the signal pulse shapes. If 
the baseband pulse duration is longer than one signaling interval, we 
assume that the symbols transmitted during different time slots are 
statistically independent. The applicability of the method to compute the 
spectral density is illustrated by examples for binary, quaternary, 
octonary, and 16-ary PSK systems with different pulse overlap. Similar 
methods yield the spectral density of the output of a nonlinear device 
whose input is a random baseband pulse train with overlapping pulses. 

I. INTRODUCTION 

In recent years, digital phase-modulation techniques have been 
playing an increasingly important role in the transmission of infor- 
mation in radio and waveguide systems. Various methods have been 
developed recently for computing spectra of a sinusoidal carrier phase 
modulated by a random baseband pulse train. 1-8 

In this paper, we derive the spectral density of a carrier phase modu- 
lated by a random baseband pulse stream in which the signaling pulse 
duration is finite and the signaling pulses may have different shapes. 
The spectral density is expressed as a compact Hermitian form in 
which the Hermitian matrix is a function of only the symbol proba- 
bility distribution and the associated column vector is a function of 
only the signal pulse shapes. If the baseband pulse duration is longer 
than one signaling interval and the pulses from different time slots 
overlap, we assume that the symbols transmitted during different 
time slots are statistically independent. 
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The present method also yields the spectral density of the output 
of a nonlinear device whose input is a similar baseband pulse train. 

The work reported here generalizes and simplifies prior results. The 
form of the present results provides an appropriate division between 
analysis and machine computation that enhances physical under- 
standing and simplifies numerical computations. We compute the 
spectra of binary, quaternary, octonary, and 16-aryPSK systems, 
with overlapping baseband modulation pulses of several shapes. 

II. M-ARY PHASE-MODULATED SIGNALS 

We seek the spectrum of the digital phase-modulated wave 

x(t) = cos[_2ttU + 0(0], fc > 0, (1) 

where 

4>{t)= t Q. k (t-kT), s h = 1,2, .j..,Jf. (2) 

A=-oo 

The discrete random process s& is assumed strictly stationary; as 
noted in (2), it takes on only integer values from 1 to M. The carrier 
frequency is f c . The signaling alphabet consists of M time functions, 
0i, 02, • • •, 9m, that may have different shapes; one of these is trans- 
mitted for each signaling interval T, to generate the digital baseband 
phase modulation <f>(t). The different signaling waveforms in (2) may 
overlap, and may be statistically dependent throughout the present 
section and Appendix A. 
For convenience, we define 9 

»(0aa«* (O ; (3) 

then 

x(t) = Re |e>' 2 "«' «(*)}■ (4) 

Appendix A shows that the spectrum of x(t) is 

p.(/) = 1P.(/ - /.) + iP.(- / " fc), 

2/ c ^|, n= 1,2,3, •••, (5) 
where 9 

P»(/) = r *.(r)r**"dr l (6) 

J — CO 

*,(t) = *„(*, r) = lim ±r \ A *,(*, r)di, (7) + 

*.(<, r) = (v(t + t)i>*(0) = (e''"^>"*W 1 ). (8) 

* The symbol denotes average on t throughout. 
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The first term of (5) is the spectrum of the complex baseband wave 
v(t) shifted to the carrier frequency +/ e ; the second term is the spec- 
trum of v*(t) shifted to — /„. The spectral relationship of (5) is strictly 
true as long as the condition on f c is satisfied, whether or not the two 
spectral terms overlap; that is, this result applies strictly to both 
narrow-band [/«. » bandwidth of P B (/)] and wideband modulated 
waves. 

The condition of (5) requires that twice the carrier frequency is not 
an integral multiple of the signaling rate l/T. P,(/) has in general 
both continuous and line spectra, the lines occurring at frequencies 
n/T for integer n. The condition in (5) guarantees that the line com- 
ponents of the two spectral terms never coincide. In the exceptional 
case, where the two sets of lines do coincide, it is not surprising that 
it no longer suffices merely to add powers of the two terms, as in (5) ; 
however, if /„ is high enough the modulated wave is narrow band, the 
two spectral terms of (5) do not overlap significantly, and (5) pro- 
vides a good approximation even if the condition is violated. 

Note that we have not found it necessary to randomize the phase 
of the unmodulated carrier or the position of the time slots of the 
digital modulation, as is done in various other studies. Thus, rather 
than (1) we might have considered 

x(t) = cos [2irf c t + <f>(t - t ) + *o], (9) f 

with <p(l) still given by (2). The spectral relation of (5) will again be 
strictly true when 2f c T = integer if <f> and t are independent, and 
either is uniformly distributed over a suitable interval (Appendix A). 
However, we retain the original formulation of (1) and (2) throughout 
— corresponding to <t> = 0, t = in (9) — to determine the effect of 
deterministic carrier and modulation phase on the statistics of the 
modulated wave. 

Consequently, we study only P„(/) throughout the remainder of 
this paper. This suffices for all cases except a low carrier frequency f c 
that is a precise multiple of half the baud rate 1/2T; the additional 
calculations necessary for this exceptional case are suggested in Ap- 
pendix A, but not carried out in detail. 

III. NOTATION AND STATISTICAL MODEL 

The introduction of vector-matrix notation greatly simplifies the 
ensuing analysis. 



* Either of the two parameters, <*>o or to, shifts the relative phase between carrier 
and modulation ; thus, only one of them is really necessary. 
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First, we rewrite the phase modulation of (2) as 



+ • • • + o^ffj,^ - fern (io) 

For a given A; (i.e., for a given time slot), one of the a k 's is unity and 
the rest are zero; 

of = 0, tV s/t, 



(U) 



where s k is the strictly stationary, discrete random process of (2), 
taking on the integer values 1, 2, • • •, M. 

Now we define for convenience the row vectors 



(12)1 



^^WaP ---an 
,g(0, = \JBi(t) 9*(t) • ■ • flfjKOl 

whose components are respectively the coefficients and pulse shapes of 
(10). The transposes 1 of these row vectors are the column vectors 



a fc ] = ,a/, = 






, g(0]=.g(0/ = 



fffiW 

ff2(0 
.0A/(O 



(13) 



Define the unit basis row vectors 



£±= [10 

e 2 = [0 1 



01 

0], 



e^s [0 ••• 1], 

with corresponding transpose column vectors ei], eaj, 
,a fc| takes on only the values gj, ^ • • • , e^ by (11) : 

Now we rewrite (10) in vector notation as 

*(fl = £ a*/g(< - M01 



(14) 

ejifj. Then 

(15) 

(16) 



where • signifies ordinary matrix multiplication throughout. The term 

t Boldface quantities denote matrices throughout. Row and column vectors are 
distinguished by the additional notation i_ i and ], respectively. 
* The transpose of a matrix is indicated by the symbol '. 
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,a.k, is a vector-valued, discrete 1 random process, strictly stationary by 
the assumed strict stationarity of s* of (2). We define the first- and 
second-order probabilities of ^or a* J as 

Wi = Pr[a* = e<) = Pr(s* = i] ^ 0. (17)* 

W„(i, j) s Pr{a* = e„ a, +n = e,} = Pr{s, = i, s k+n = j\ ^ 0. (18) J 

Wi is the probability that the ith signaling waveform g t is transmitted 
in any time slot, W „{i, j) the joint probability that the signaling 
waveforms 0,- and g, are transmitted in two time slots separated by n 
signaling intervals. Wi and W„(i, j) are independent of k by the assump- 
tion of stationarity. Then since the marginal probabilities are obtained 
by summing over the joint probability function, 



M M 

E W»(i, j) = Wj, E W n (i, j) = Wi. 

i=i i=i 

Normalization of the total probability to 1 requires 



E Wi - 1, 



at m 

E EF.ft;] = l. 
i=i y=i 



(19) 



(20) 



Now we introduce vector-matrix notation for the probabilities. Let 

,£,= [wi Wi • ■ ■ w M ~\ (21) 

be the probability row vector whose elements give the probabilities of 
the different signaling waveforms, with transpose column vector 

w]= l w J '. (22) 



Let 



W„ 



Wn{l, 1) 

^n(2, 1) 



W n (l, 2) 

W n (2, 2) 



W n (l, M) 

W n (2, M) 



(23) 



W n (M,M) 
oint probabilities of all 



W n {M,\) W n {M,2) 

be the matrix whose elements specify the 

pairs of signaling waveforms separated by n time slots. Further, define 

J, - 1 J ■ D 1 • ■ • 1] (24) 

as a vector with all M elements unity. Then (19) and (20) may be 

J -W n = l w J , W n -l] = w] (25) 



written as 



t.ajt.is defined only for integer values of its independent variable k. 
1 It is understood that a, e are either row or column vectors throughout (17) and 
(18), and in succeeding equations. 
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and 

JL . w ] =i w i .l]= 1, !-W n -l] = l. (26) 

The probability matrices have the following useful properties. For 
n = 0, the joint probability matrix of (23) and (18) is diagonal, with 
diagonal elements the first-order symbol probabilities; 





Wo = 



Wl 



Wi 



= Wd, 



(27) 



w M _ 

where we define the diagonal matrix w d for later convenience. Then 
Wrf-1] = w], i-Wd = ffl. (28) 

By symmetry, 

W n (i, j) - W- n U, i), (29) 

or in matrix notation 

W n = W_ n '. (30) 

We assume that signaling waveforms in widely separated time slots 
are statistically independent: 

(3D 



lim W»(i, j) = Wi Wj, 

n-»» 



lim W n = w]-ffl. 



(w]-ffi)' = w]-a, 



(32) 



(33) 



or in matrix form 

Since 

(30) and (32) yield 

W„ = W„' = W_„. (34) 

Using the above, the mean and covariance of the strictly stationary, 
vector-valued, discrete random process a* are 

(a*]) = w], (&> = ffi> (35) 

O a (n) = (**+•"}'&) = W„. (36) 

The assumption of independence in (31) and (32) imphes uncorrelation 
in (36) as n — * °o . Conversely, uncorrelation implies independence, 
because the random vectors a* are unit basis vectors (15), 10 rendering 
the covariance and probability matrices identical. Thus, rather than 
assuming independence in (31) and (32), we might equally well 
assume that the modulation vectors & k become uncorrelated for 
widely separated time slots. 
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The spectral density of a* is the discrete transform of (30) ; 

P.(/) - £ e-*"»W„. (37)t 

n= — » 

Equations (36) and (37) are the matrix extensions for discrete vector 
random processes of similar scalar relations for discrete scalar random 
processes; 11 the symbol ~ indicates that a discrete random process is 
under consideration. We separate (37) into line and continuous spectral 
components : 

P.(/) = Pa,(/) + Pa c (/)- (38) 

Then from (32) and (37) : 

P.,(/) = w]- a £ d(f-n), (39) 

P..(/)= £*-*""{ W.-w]-ffi}- (40) 

n=— co 

The assumption of (31) and (32) that signaling waveforms in widely 
separated time slots are independent, or the equivalent assumption 
that a*+„ and at become uncorrected for large n, eliminates line 
components except dc if we confine our attention to the fundamental 
frequency interval |/| < 5, and so retain only the n = term in 
(39). * Consequently, the modulation has no periodic patterns. 

IV. PSK AS A BASEBAND PULSE TRAIN 

Our problem has been reduced to determining the spectral density 

of the complex wave v(t) : 

v{t) = e*<», (41) 

<t>(t) = £ £h>-g(t ~ mi, (42) 

k =-00 

where a and g are M-dimensional vectors. We show that if the signal- 
ing pulses are strictly time-limited to an interval KT, v(t) may be 
written 

v(t) = t fafr(« - kT)l. (43) 

A =—00 



+ While this relation is defined for the entire range - » < / < °° , Pa (/) is periodic 
in / with unit period, and only the fundamental period |/ 1 < | is normally of interest. 
Thus, the inverse transform is 

*.(n) = W„ = J* Pa(/)e+»""d/. 

* We may thus write a* = w + a c i with (40) giving the spectral density of a c and 
(39) the spectral density of w. 
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Vectors b and r are M x -dimensional vectors expressed in terms of a 
and g, respectively. The different terms in (43) are strictly nonover- 
lapping: 

r(fl] = 0], t£ 0, t> T. (44)* 

Each b k is a unit basis vector, i.e., only one of its M K components is 
unity, all others being zero. The spectral density of (43) is determined in 
Appendix B. 

Figure 1 shows portions of <l>(t) of (42) for four different maximum 
signal pulse durations; the terms k = - 1, 0, 1, 2 of (42) are shown, 
and for convenience a k has been taken the same for each of these time 
slots. The pulses are positioned along the time slots such that the 
limits of each signal pulse lie on the boundary between adjacent time 
slots (i.e., t = integer- T); since symmetric pulses have been chosen for 
illustration in Fig. 1, their maxima are centered in the time slots for 
K odd, and lie on the time-slot boundaries for K even. Examine the 
(0, T~\ time slot in Fig. 1 as typical; then with the above choice of 
pulse positions, the number of pulses contributing to 4>{t) in each 
time slot equals K. Since each pulse can take on M different shapes, 
<t>(t) can take on M K different shapes in each time slot; the same is 
true for v(t) of (41), thus demonstrating the representation of (43). 

It remains for us to express the pulse shapes r(t) and coefficients 
bfc of (43) in terms of the signal pulses g(0 and coefficients a. k of (42). 
We give separate treatments for the cases K = 1 and K = 2, and 
extend these results to general K. 

4.1 Nonoverlapping pulses: K = 1 

The top portion of Fig. 1 shows digital phase modulation for which 
the signal pulses in different time slots never overlap; in this case, 

g(0J -01 t£ 0, t > T, (45) 

where 0] represents the M-dimensional zero vector. Define 

jVftCt) e ta(D . . . ei».w(')] j o < t^ T. 

q(0 - (46) 

& t£ 0, t> T. 

Then (41) and (42) may be written 

v(t) = t s*,-q(< - M91 ( 4? ) 



t o] = l J ' is a vector of appropriate dimension (here M K ) with all elements zero. 
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K = 4 



Fig. 1 — Phase modulation for different signal pulse durations. Index k [eq. (42)] 
is shown near peak of each pulse. Also, for simplicity, same signal pulse is shown for 
each k. T = time slot duration or signaling period. KT = maximum signal pulse 
duration. Note different pulse center location for odd and even K. 



Comparing (47) and (43), the parameters of the latter are given as 
follows for nonoverlapping signal pulses: 



>*, = £*, r(0] = q(0]; K= 1. 



(48) 



4.2 Overlapping pulses: K = 2 

This case is illustrated in the second portion of Fig. 1. In the (0, !F] 
time slot, the k = 0, 1 pulses contribute. We have 

f(fl] = 0], t £ - T, t> T. (49) 
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Define 

[ e fo<*> eiatW . . . e *«(D] } —T < t ^ 7\ 

q(0 - (50) f 

& t^ - T, t> T. 

Then 

w(0 = £ {^-qtf - /cT)]} fcu+i.-qft - (ft + l)? 7 )]} 

= £ f.a*.-q(* ~ kT) ]} X |.a fc+v q(t - (ft + DO 

ft— 00 

- £ l&X fei Mq(<-fc?')]Xq(/-(fc + l)?')]), (51) 

* =— 00 

where X denotes the (right) Kronecker product 1213 throughout. 1 
Comparing the last line of (51) with (43), the parameters of the latter 
are given as follows when no more than two signal pulses overlap: 

b*, = 1 a tl X I a i±v r(<)] - q(fl] X q(( - D]; K = 2. (52) 

Since the elements of b* consist of all pairs of products of the elements 
of a fc and a*+i, and since the a* are 2l/-dimensional unit-basis vectors 



+ Comparing (50) with (46), note that the definition of q(0 is different for different 
K; q(t) ^0 over the same interval in which g(J) may be nonzero. 

* The second line of (51) follows from the first by the observation that the two 
scalars may be regarded as 1-by-l matrices ; consequently, their scalar product and 
their Kronecker product are identical. The third line follows from the second by the 
well-known result connecting ordinary matrix and Kronecker products as follows: 
Consider two arbitrary matrices A and B with elements a,-,-, by. Define the (right) 
Kronecker product by (Refs. 12 and 13) : 



AXBe 



[eiiiB ai 2 B • • • "I 
a 2 iB a 22 B 



The following results may be found in Refs. 12 or 13. 
For any matrices A, B, C, and D : 

A X B X C = (A X B) X C = A X (B X C), (A X B)' = A' X B'. 

For A and B the same size, and C and D the same size (possibly different from the 
size of A and B) : 

(A + B)X(C + D)=AXC + AXD+BXC + BXD. 

Assume the matrices Ai, Bi and A 2 , B 2 are dimensioned so that the ordinary matrix 
products At-Bi and A 2 B 2 exist (i.e., there are \ columns of Ai and rows of Bi, and 
M columns of A 2 and rows of B 2 ). Then 

(Ai-Bi) X (A 2 -B 2 ) = (A, X A,)-(Bi X B,) ; 
this result generalizes to 

(Ai-Bi-C) X (A 2 -B 2 -C 2 ) X (A 3 -B 3 -C 3 ) 

= (A, X A 2 X A 3 ) • (Bi X B 2 X B,) • (Ci X C 2 X C) 
etc. 
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by (14) and (15), the b* arc A/^-dimensionaJ unit-basis vectors: that is, 
b k has one element unity and the remainder M 2 — 1 elements zero. 
Note from (52) and (50) that 



r(0] = 0], t £0, t> T. 



(53) 



We illustrate these relations for the binary case, 14 in which only 
two signal waveforms gi(t) and </ 2 (0 are transmitted. Then 



-TO-ISH!!]' ^-T, t >T. 



q(0] = 



[2] 



- T < t ^ T. 



r(0] - 



e Jla,U)+ 0l (t-T)] 
e ilg,U)+0iU-T)] 
e j{ui (')+<;, U-T)] 
e Jla 2 (i)+atU-T)] 



t£ - T, t> T. 



< t ^ T. 



t ^ 0, t > T. 



h, = w op: x cofti off.,: 

= [a* 1 * aft, <*> afti «P «fti a * 2) "fliJ 





[1 0] 


[0 1] 


[1 0] 


[10 0] 


[0 10 0] 


[0 1] 


[0 10] 


[0 1] 



(54) 



(55) 



(50) 



(57a) 



(57b) 



The four possible waveshapes for <f>(t) in the (0, T~\ time slot are 
shown in Fig. 2. 

4.3 Overlapping pulses: General K 

The general case is treated by straightforward extension of the above 
method; Fig. 1 illustrates the phase modulation for K = 3, 4. The 
following expressions differ .slightly for the even and odd cases; they 
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9l<t) 




g,(t-T) 

\ 

\ 
\ 

\ 


91 (t) 

/ 
/ 
/ 


Y 

\ \ 

\ ^ 




s 

/ 
/ 
/ 

y 


X 









T 
g,(t-T) 





\\ 


/ 
/ 

g 2 (t-T) 





// 
// 
// 


\ 
\ 
\ 

T V 







T 


\ 

\ 
\ 


/ 
/ 
// 




\ 

\ 
\ 
\ 
\ 

92 (t) 


X 


/ 
/ 
/ 


92 (t) 




g 2 (t-T) 



Fig. 2—^(0 for < t ^ T, binary signaling, and K = 2. 



correctly reduce to the above results for K = 1, 2. 

g(o: = o], 



t£ ~^±T, *>^±ir; if odd. 



(58) 



li-fr, t>^T; 



K even. 



K- 1 



T <t^ ^4^ Ti K odd. 



-fr<«*fr. 



if even. 



&(*), - 



(59) 



t ^ - ^^ r, t > ^^ T; K odd. 



A 



'*-!*, «>f r ' 



if even. 
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The parameters of (43) are: 



(K-D/2 (K-D/2 

l bfc J = IIx to '(<)]- IIx q(«-*TO; ^odd. (60)* 

i=-(A-l)/2 t=-(X-l)/2 

^= ffx fi*±i, r (0]= IIx q(*"<r)]; #even. (61)* 

.=-(A72)+l i=-(iC/2)+l 



Note that 



r(0] = 0], tgO, t>T. 



(62) 



V. PSK WITH NONOVERLAPPING, CORRELATED SIGNAL PULSES: K = 1 

Assume the signal pulses are confined to one time slot, as in (45) : 

g(0: = 0], t£ 0, <>r. (63) 

From (46) to (48) and (160) to (161), 



R 



(/)]= r 

Jo 



S-JIr/l 






tft. 



(64) 



We separate w(<) of (2) and (3) or (41) and (42) into line and continuous 
components: 

(Q = e i*(o = W| (i) + Vc (t). (65) 

Using (48) and comparing (32) to (37) with (149) to (153), we have 
from (171) 



J 



(«] 



From (165), 



vi(t) = yffl. Z R( 5]|a*-*". 



>..(/) - f 2 lffi-R(/)]l 2 n f m 5 (/~|)' 



(66) 



(67) 



* The symbol IIx used here and subsequently indicates a multiple Kronecker 
product : 

N 

IIx Ai ■ Al X Al+i X X A w _, X A*. 

i-L 

The products in (60) and (61) contain K factors. 
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which may be written 

P.,(/) = ^ \wiRi(f) + »««»(/) + ••• + w M R M (f)\ 2 

From (162) and (164) 

P..(/) = 4 .R(/), - £ ^ Jlr/r " (W„ - w]-wj -R*(/)]. (69) 

i N = — CO 

We illustrate these general results with two examples. 

5.7 Independent, M-ary signal pulses 

Equation (32) now holds for all n ?* 0; using (21), 



W„ = w]-ffis 



t0jf 



[i<;i iu 2 • • • wm~], n^O. 



(70) 



From (27), 



Wo = w„ = 



w j 



10a 



Wtf 



(71) 



Then only the n = term remains in (69) ; 

P..(/) = ^RCQ-w«i-R*(/)] - ^ [s-R(/)]| 2 . 

Writing out the matrix operations, (72) yields 

P..(/) =^[^il-Ri(/)l 2 + ^|fl 2 (/)| 2 + ••• +^M|«Af(/)| 2 ] 
- i |i0ii2i(/) + wMf) + ■ • • + Wjrfiji(/)| a 

1 M M 



(72) 



(73) 



This result has been given in Ref. 6. The line component is given by 
(66), its spectrum by (67) and (68). 
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5.2 Correlated, binary signal pulses: M — 2 

The matrix W n consists of the four upper left-hand elements of 
(23). The constraints of (19) and (20) or (25) and (26) render three of 
these functions dependent on the fourth. 

W n = [W n (\, 1) - ioX] \\ ~J] + w]. ffi . (74) 

From (27) 

Wo(l, 1) = wi, (75) 

from (30) 

W n (l, 1) = W- n (l, 1), (76) 

and from (31) 

lim W n (l, 1) = wl (77) 

n-*oo 

Equation (69) becomes 

P*(/)- h,\Ri(f)~RM\ 2 t e-*"T"[W n {l,l)-w\-]. (78) 

/ n= — oo 

W n (l, 1) may be any discrete covariance function satisfying the con- 
straints of (75) to (77). The line component and its spectrum are again 
given by (66), (67), and (68) with two-component vectors, e.g., 



P.i(/) = ?-* \wiRi(f) + mR*(f) I 2 n t m *(/'|)- 



(79) 



The binary independent case is obtained by setting 

W n (l, 1) = w\, n9*0. (80) 

Then (78) becomes 

Pv e (/) = ^- 2 |fli(/)~ R*(f)\ 2 - (81) 

This agrees with (73) for the binary case. 

VI. PSK WITH INDEPENDENT, OVERLAPPING SIGNAL PULSES: K = 2 

Assume that no more than two signal pulses overlap, as in (49) : 

g(fl] = 0], t^ - T, t> T. (82) 
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R(/) is determined from (50) to (53) and (160) to (161): 



R(/)]= f 

Jo 



»-j'2ir/* 



eiloiM+oud-T)) 

e ilo 2 (.t)+g u (.t~T)} 
e i{o 3 (t)+ 0l (.t-T)) 

e J{0M-i(*)+0u<.t-T)) 

e i[oH(t)+ 0l <,t-T)\ 



dt. 



(83) 



From (52), 



(84) 



We determine the mean and covariance of b* from (32) to (37) and 
(149) to (153). We assume throughout this section that signal pulses 
in different time slots are independent, as in (70) and (71) : 



Wo = w d ; W n = w]-ffi, \n\ ^ 0. 
For the mean, since a* and a.k+i are independent, 

= ,w' 2 J. 



(85) 



(86) 



We introduce the notation that an integer exponent enclosed in square 
brackets denotes the Kronecker power, 12 i.e., the Kronecker product 
of the matrix (or vector) with itself the indicated number of times. 
For the covariance (see footnote, p. 908) : 



Since 



Ob(w) m (bt+J-^fc,) 

= <(a*+n] X a fc+ „ + i]) • (a*, X ,a*+0 ) 
= ((a-k+nl-Qk) X (a fe+n+ i]-,a fc -n)). 

5» b (n) = & b (-n), 



(87) 
(88) 



we study only n ^ 0. We treat three cases. 

First assume \n\ ^ 2. All the a's in (87) are independent; from the 
second line, 

* b (n) = (w]X w])-(fflXffi) 

= w]' 2 '-^ 2 ', \n\ ^2. (89) 
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Next consider n = 0. From the third line of (87) and the fact that 
a* and at+i are independent, 

Ob(0) = wS 2] . (90) 

Finally, for \n\ = 1, from the third line of (87) we have t (see 
footnote, p. 908): 

* b (l) = <(a* + i] X a^) X (a t+2 ]X.a t+ Q) 

- <&, X (a* +1 ] X a^+i) X a* +2 ]>. (91) 

Since a*, a*+i, a fc+2 are independent, 

*b(D = a X w d X w]. .(92) 

From (88), 

* b (-l) = w]X w d Xw;. (93) 

From (8G) to (89), (153) is satisfied, i.e., D* for widely separated 
time slots are uncorrected. * Therefore, from (86), (165), and (171): 

vi(t) = f& [2] - n £ a R (%)] ein2r " T '> < 94 ) 

P.«(/) - 4lffl m -R(/)]l 2 t KS-nT). (95) 

1 n = — oo 

In comparing these results with (66) and (67), note that R(/) of (83) 
differs from R(/) of (64). 

Substituting (86), (89), (90), (92), and (93) into (162) and (164), 
the continuous spectrum is 

+ (ffi X w d X w] - w2 l2] -SL w )e- i2T/T 
+ (w] X w d X ffi - w]™ -&w)e+i 2 " T } -R*(/)]. (96) 

We illustrate these results for the binary case. 



* The following vector relations may be established by inspection: 

a]-«b,= a] Xfe-JbXa} 

1 bk+n and b* also become independent as n — * », because the b* are unit basis 
vectors (Ref. 10). 
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6.1 Independent, binary signal pulses: M = 1 

R(/) in (83) now contains four components, the Fourier transforms 
of the components of (56). We have 

^[2] = [yfi Wl W 2 WlW 2 W2~\ 





w\ 


0' 




wS a] = 


w\W 2 
W\Wz 




*-. 


wl 






f w\ w\w 2 w\w 2 w\wV\ 


W ]I2]. ffi [2] = 


w\v)2 wiwl w\w\ W\w\ 

W\V)2 W\w\ V)\w\ W\U)\ 




wiwl W1W2 iviwl w\ 




w\ wiw 2 


ffi X w d X w] = 


wiiv 2 Wiwl 
wiw 2 wiu 






Wiwl 


wl 



w] x w d x a = fex^x w])'. 



(97) 
(98) 



(99) 



(100) 



(101) 



(102)1 



Equations (98) to (102) substituted in (94) to (96) yield the final 
results for independent, binary signal pulses. 

To illustrate, we write out a few terms of the matrix contained in 
{ } in (96), for the continuous spectrum in the binary case: 15 



}n = wi{\ - wi){l + 2wi cos 2irfT), 

} H = wiw 2 e + **f T - wiw 2 (l + 2 cos 2tJT), 



(103) 



Writing out even the present simplest overlapping case produces an 
awful mess. Consequently, in the examples presented in this paper the 
digital computer is programmed to work directly with (96) or its 
generalization (116), performing matrix operations (both ordinary and 
Kronecker matrix multiplication) directly, rather than entering expres- 
sions such as (103). In this way, quite complicated cases involving 
multilevel signal pulses overlapping several time slots may be simply 
treated. 



See footnote, p. 908. 
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The general results for K = 2, (94) to (96), require the signal pulses 
to be less than two time slots in duration (82). Therefore, they must 
apply when the signal pulses are restricted to a single time slot, i.e., 
when the stronger condition (63) is satisfied, and must then reduce to 
the results for K = 1, (66), (67), and (72). This reduction is demon- 
strated in Appendix C. 

VII. PSK WITH INDEPENDENT OVERLAPPING SIGNAL PULSES: GENERAL K 

Finally, assume that the signal pulses occupy at most K time slots, 
so no more than A' signal pulses overlap ; g(t) is now given by (58). 
r(0 is given by (59) to (62), and R(/) by (160) and (161). 

The b k are given by (60) or (61); the mean and covariance of b* 
are found from (32) to (37) and (149) to (153), assuming throughout 
that signal pulses in different time slots are independent, as in (70) 
to (71). 

For the mean, 

,6 = (b*> - IIx (**+<) (104) 



by the independence of the a*. The limits over the II x> given in (60) 
or (61) for K odd or even respectively, extended over K factors in either 

case. Thus, 

3 = wi*i. (105) 

£j i — i 

For the covariance (see footnote, p. 908) : 

* b (>0 - (b*+n]-J>&) 

- <(IIx **+»+.-])• (II x2*±i)> 

= aix(a* + » + .-]-a^)>, (106) 

i 

the IIx being taken over A' factors as given by either (60) or (61). 
By stationarity, (106) is independent of k; consequently, for both even 
and odd cases 

O b (n) = (ba+iJ'lu) 

= ((nxa„ + ,])(nx^)) 

= (nx(a» + .]^)- (107) 
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First, for n = 0, from the third line of (107) and the independence 
of the at, 

*b(0) = fix <»<]-.§*> = w ^ 3 - (108) 

Next consider n = 1. From the third line of (107) (see footnote, 
p. 915), 



*b(l) = (nx(faXa m ])\ 



= fei> X j n x <a,]-^)| X <a* +1 ]> 

= ffi X wS x - 1] X w], (109) 



the second line again following from the independence of the a*. 
For n = 2, proceeding as above: 

*b(2) = /fixCaiXa^])^ 

-= <&]) X ^ftx (a.- + i] X^ X (ajc + .]> 

- <&l> x ^ fix <&, x a <+iH)) x <■**■]) 

= <ai) X <aj,) X 



II x (a.H'&M X (ajc+J) X (ajc+J) 

=,w/ 2 i X wS M] X w]™. (110) 

By induction: 



*b(3) = ffi f3] X wS K - 3] X w]") 



(111) 



*„(£ - 1) = a'*- 1 ) X w d x w]'*- 1 ] 
•kW = ffi m X w]m. 
Next, from the second line of (107) and the independence of the &k, 
*b(n) = WJW-ffi'fl, n ^ if. (112) 

This result is of course consistent with the final line of (111) (see foot- 
note, p. 915). Finally, 

*b(n) = &(-n). (113) 
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From these results the line component is : 

vi(t) = Y &lK] \j^„ R ( r f)] ein2,rt ' T; (114) 

P»,(/) = j- 2 lffi [K1 -R(/)]l 2 n f M 8 ( f ~l')' (115) 

The continuous spectrum is : 

P.,(/) = f.R(/H(*S* ] - wjm.E 1 * 1 ) 

+ E* (ffi 1 " 1 X wS K ~" ] X w]i"i - w]f*J-E(*])e- i ' ,2 ' r/7 ' 



K-l 

E 

n = l 
K-l 

E 

71=1 



+ E (w [n] x wS K - n] x a 1 " 1 - w2 lK] -sL lK} )e^ inl * fT ] 

•R*(/)]- (H6) 

If we set K = 2 in (114) to (116), these results agree with those of 
Section VI. If the signal pulses are restricted to a smaller number of 
time slots K < K, the present results reduce correctly to those ap- 
propriate for K; this reduction is shown for S. = 1, K = 2 in Appendix 
C, but the general case is left as an exercise for the reader. 

VIII. EXAMPLES 

We now consider a number of examples of computing spectra of 
multiphase PSK systems, subject to the following assumptions: 

(i) The number of phase levels is a power of 2, 

M = 2 N , N an integer. (117) 

(it) The M signaling pulses have a common shape, and the M 
phase levels (peak phase modulations) are equally spaced. 

iS-jD 3 • • • (2M - l)]ff(0, (H8) 

<7(*) m a* = I- (H9) 

Figure 3 shows vector diagrams of the peak phase modulations 
for M = 2, 4, and 8. We also assume that g(t) is symmetric, 

g(t) = g(- t), if even, 

g(t + T/2) = g{-t+ T/2), K odd, ^ ' 
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m = 2 




m = 8 

Fig. 3 — Peak phase modulation for PSK with equally spaced levels. 

with maximum at 

t = 0, K even, 

t = j, K odd. 



(121) 



(Hi) Signal pulses in different time slots are statistically indepen- 
dent, and all signal pulses are equally likely; 



w] = ^ 1], w d = ± I, 



(122) 



where I is the identity matrix of order M. 

As a consequence of (118) and (122), we can show that P„(/) is 
symmetric, that is 

P.(/) = P.(-/). (123) 

8.1 Rectangular nonoverlapping signal pulses 

If g(t) is a rectangular pulse (see Fig. 4) of duration r? ^ T, 



0, otherwise. 



(124) 
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From (118) and (64), 



Ri(f) = 



expjj 



M 



&-!)}- A s ™^fi><ri*fT 



irf 

TTJ 



M. (125) 



^ , < t^ T 



The spectral density P„(/) can be computed from (68), (73), and 
(122) with R,(/) given by (125). The results can be shown to agree 
with those given in Refs. 4 and 6. Since this case can be treated as 
amplitude modulation and has been discussed extensively in Refs. 4 
and 6, we shall not discuss it further in this paper. 

8.2 Raised-cosine nonoverlapping signal pulses: K —1 

If g(t) is a raised-cosine pulse (see Fig. 5) of duration T, the pulse 
just fills the time slot, and 

IT 2irt~\ 

^l-cos^-j, 

g(t) = (126) 

0, otherwise. 

In this case, R(/)] may be evaluated either numerically or using 
Bessel function expansion (which again requires the use of a computer). 
For M = 2, 4, 8, and 16, we have evaluated as above 

P.U) = P,(/) + P«m (127) 

with results shown in Fig. 6. We note that there is very little variation 
in either the discrete or the continuous spectrum when M is increased 
from 2 to 16. The tails of the spectrum are not shown in Fig. 6, although 
they are easily calculated down to the — 60-dB level. 



9klt) 






T-T) 

2 



1^1 
2 



Fig. 4 — Square-wave signaling of duration rj g T. K = 1. 
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Fig. 5 — Raised-cosine signaling with pulse duration T. K = 1. 

8.3 Raised-cosine overlapping signal pulses: K = 2 

If a raised-cosine signal pulse (see Fig. 7) just fills up two time slots, 

irt 



|ri + cos|l, -T<t^T, 



g(t) = (128) 

0, otherwise, 

K = 2, and two signal pulses overlap. 

In this case, P„,(/) and P„ e (/) are given by (95), (96), (122). Using 
a digital computer, P„(/) has been determined for M = 2, 4, 8, and 
16, and the results plotted in Fig. 8. Again, it may be noted that there 
is not very much variation in either the discrete or the continuous 
spectra when the number of phase levels M is increased from 2 to 16. 

Comparing Figs. 6 and 8, we observe that the power in the line 
components in a PSK system with overlapping signal pulses is smaller 
than the power in the lines with nonoverlapping pulses. Also, for the 
same signaling rates, the principal portion of the continuous part of 
the spectrum with K = 2 is narrower than that with K — 1 . The tails 
of P„(/) were calculated down to —60 dB even though they are not 
shown in Fig. 8. 

In Fig. 9, for M = 4, and raised-cosine wave signaling, we plot the 
spectral density P»(/) when the amount of overlap is increased from 
zero to one time slot. Specifically, we assume that 

I Tl + 608^1, -(3T<t^(3T, 0.5^/3 5; 1, 

g(t) = (129) 

0, otherwise. 

Note that K = 1 when = 0.5, and K = 2 when 0.5 < ^ 1. 
Figure 9 shows that there is a gradual reduction of power in the line 
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Fig. 6 — Spectral density of binary, quaternary, octonary, and 16-ary PSK systems 
with raised-cosine signaling and pulse duration T. K = 1. Although the values of the 
spectral density for M = 2, 4, 8, and 16 are shown to be the same over certain seg- 
ments of the range of /, they are usually different from each other, but this difference 
is not large enough to be shown on the linear scale in Fig. 6. 



components when the overlap is increased from zero to one time slot. 
There is also the narrowing of the principal portion of the spectrum. 

Decrease of the carrier component [n = term in (08) or (95)] when 
jS is increased may result from the fact that the phasors in Fig. 3 spend 
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Fig. 7 — Raised-cosine signaling with pulse duration 2T. K = 2. 

increasingly less time in the neighborhood of 0°. There are also smoother 
transitions between the phasors with increasing 0; this probably ex- 
plains the decrease in width of the main part of continuous spectrum. 
The narrowing of the spectral density and the reduction of the discrete 
spectral lines may or may not indicate better interchannel and inter- 
symbol performance, but this remains a subject for future investigation. 

IX. SUMMARY AND CONCLUSIONS 

Matrix methods have been used to express the spectral density of a 
sinusoidal carrier phase modulated by a pulse train with a finite pulse 
duration. Arbitrary pulse shapes may be used for M-ary digital signal- 
ing and they may overlap over a finite number of signaling intervals. 

If the pulse duration is one signaling interval, our results give the 
spectral density even though successive symbols are correlated. If 
the pulse duration is more than one signaling interval, we assume that 
symbols transmitted during different time slots are statistically inde- 
pendent; the case of correlated symbols may be considered by ex- 
tension of the present work, but the required statistical description of 
the modulation a k will be more complicated. For example, if K = 2, 
in addition to (a*]), (a*] X &i), we need to know <a,] Xa t ]X^X a™), 
the fourth order statistics of a*, to determine the spectral density from 
our methods. 

The spectral density has a compact Hermitian form suitable for 
numerical computation by a digital computer. The associated Hermi- 
tian matrix is a function of only the symbol probabilities, and the 
column matrix associated with the Hermitian form is the Fourier 
transform of certain time functions related to the signaling pulses. The 
computations presented in this paper for binary, quaternary, octonary, 
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Fig. 8 — Spectral density of binary, quaternary, octonary, and lG-ary PSK systems 
with raised-cosine signaling and pulse duration IT. K = 2. Although the values of 
the spectral density for M = 2, 4, 8, and 16 are shown to be the same over certain 
segments of the range of /, they are usually different from each other, but this dif- 
ference is not large enough to be shown on the linear scale in Fig. 8. 
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Fig. 9 — Spectral density of quaternary PSK system with raised-cosine signaling 
and pulse duration 2/32', 0.5 ^ ^ 1. When = 0.5, K = 1, the pulse fills up just 
one time slot and the spectral density is as shown in Fig. 6. When = 1, K = 2, 
the pulse fills up just two time slots, and the spectral density is as shown in Fig. 8. 
Although the values of the spectral density for = 0.5, 0.571, 0.667, 0.8, and 1.0 
are shown to be the same over certain segments of the range of /, they are usually 
different from each other, but this difference is not large enough to be shown on the 
linear scale in Fig. 9. Note that when = 0.5, 0.571, 0.667, 0.8, and 1.0, 1/0 = 2, 
1.75, 1.5, 1.25, and 1.0. 
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and 16-ary PSK systems do not consume very much computer time 
and are quite inexpensive. 

Extraction of a timing wave is often essential in the detection and 
regeneration of digital signals. In a self-timed digital system not 
containing a timing wave, the wave is extracted from an incoming 
pulse stream by a nonlinear processing of the signal. For example, the 
incoming pulse stream may be passed through a square-law rectifier 
and a linear narrow-band filter to get the timing wave. We would like 
to note here that the methods given in this paper can be extended to 
determine the spectral density of the output of a nonlinear device 
whose input is a random time pulse train of the form (16) or (42). The 
results presented here apply to the particular nonlinearity exp j(-); 
other nonlinear functions are treated in a similar manner. 
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APPENDIX A 

Spectra of Complex and Real PSK Waves 

From (1) written in complex form, the covariance of the real wave 
x{t) is 



*»(r) = *,(«, r) 



= }{0J»*/«rft a ( r )-l-4r**/«'$;( r ) + e'' 2 ''« r e*''«'$„»(«, r) 



+r-*»*'«r*»/.*«i.(|, *)}, (130) 
where the cross-variance is given by 

*■••(*, r) = (v(t + t)w(O) = (•'!♦««•*>+♦<«">. (131) 

Now <I> „„•(£, t), with t regarded as a parameter, is periodic in t with 
period T; 

$„*(/ +T,t) = *„.(*, r). (132) 

This follows from the assumed strict stationarity of the s* of (2) . Then 
we may write 

*..•«, t) = £ Vn (r)e^ M ' T . (133) 
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Substituting in the time average quantity of (130) and interchanging 
the order of summation and integration, 

«*"•'*...«, r) = t <p n (r) lim ± f A a**Wrf-C/ni« 

- , . ,. sin27r[2/ e + (n/r)]A 

= .£. ^ n(T) 1™, 2x[2/ c + (n/D3A (134) 

The lim — ♦ if 2/ f + (n/T)^ for every integer n, i.e., if twice the 

carrier frequency is not a precise multiple of the modulation fre- 
quency 1/T. Under this condition, (130) and (134) yield (5). 
Next, assume that P,(/) is strictly band-limited; 

P.(/) = 0, |/| fc /,. (135) 

Then P,(/ - / e ) and P,(-/ - /„), the two terms appearing in (5), 
do not overlap. Now define 

v (t) = a*"-«i>(t). (136) 

Then 

<*>,(*, t) = e*"« r *,(*, r); *,(r) = e' 2ir '"«S> B (r). 

P.(/) = P.(/ - /c); P^O) - P.(-/ - A)- 

Therefore, subject to (135) 

/£ 0. 

P> (/) = o, 

/ ^ 2/c. 

/ ^ - 2/ c . 

P>'(/) = o, 

/£ 0. 

Now the cross-variance of j*(i) and v*(t) is 



(137) 



(138) 



$„,( T ) = $„♦(*, r ) = «i»-/.r «**/•«*„.(*, r), (139) 

where $ B „* is given by (131). The Fourier transform of (139) is the 
corresponding cross-spectrum P„,*(/). Since the two spectra P,(/) and 
P,*(/) do not overlap by (138), the cross-spectrum P w *(/) must be 
identically zero, 16 and hence also the cross-variance ; 

# w .(r) - 0, P,(/) =0 for |/| ^ /„. (140) 
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Substituting (139) and (140) into (130) and taking the Fourier trans- 
form, 

P.(/3 = IP.(/ - fc) + 1P.(-/- /«), 

P.(/) = for | /| ^ / e . (141) 

While P„(/) will never be strictly zero, it will eventually fall off so 
rapidly with increasing |/| that the spectral relation of (5) will be a 
good approximation when f c is large enough so that the two terms do 
not overlap appreciably, even if 2/ c is an integral multiple of l/T. 

Consider now the exceptional case of a low carrier frequency that is 
an integral multiple of half the modulation frequency: 

2/ c + ~° = 0. (142) 

Then in (134) the lim — > 1 for the n term, and -* for all other 
terms as before. Therefore, (130) and (134) yield 

*,(r) = ie^'O/r) + <p-2/ c t(t)1 

+ \e-v< T L$l(T) + vU fcT (r)], 2f c T = integer. (143) 

The exceptional case requires the evaluation of the additional quantity 
<P- 2 f eT (T), 2f e T = integer, where <p is denned in (133). This may be 
done similarly to the evaluation of $„(t) performed in the text, but 
it is not undertaken here. 

Finally, consider the wave of (9), with <f>(t) given by (2), with the 
carrier and modulation phases </> and k independent random variables, 
independent of the modulation. Equation (130) is replaced by 

<*>*(r) = i{e> 2 *'"* B (r) + e-^f^Kr) 



+ e **fcT( e J24> a )( e J*Tj c l ) e J4*/ct$ vv ,( ti T ) 



+ e -K*fcT( e -X4 a )( e -ii*fcu) e -ji*fct^(t ) T )] } (144) 

where v, $„ and *„• remain as given in (3), (7), (8), and (131) for 
t = 0. The last two lines of (144) vanish if any of the three quantities 

«*"«'*„.(<, t), (e**°), <e*"«'.) (145) 

vanish. Therefore, 

P.(/) - iP.(/ - /«) + ».<-/ ~ /.) ( 146 ) 

DIGITAL PM SPECTRA 929 



if any of the following conditions are satisfied: 

2f c T ^ integer. 

2f c T = integer, <f> is uniformly distributed over an interval 

it, 2r, 3tt, •••. (147) 

2f c T = integer, t is uniformly distributed over an interval 

T T T 

2 tt4^ , 3 



2f c T' "2f c T' J 2f c T' ' 

The first condition is, of course, that of (5). The last two conditions 
show that suitably randomizing the phase of either the carrier or the 
modulation suffices to yield the simple spectral result of (5) when 
2f c T = integer. 

APPENDIX B 

Spectrum of a Baseband Pulse Train With Different Pulse Shapes 

We determine the spectral density of (43), 

v(t) = £ b^T(i - kT)l (148) 

by the vector analog of the corresponding derivation for a train of 
equally spaced pulses with similar shapes. In this appendix, the wave- 
forms r for different time slots (i.e., different k) may overlap, although 
they do not overlap in the applications of this paper. Also, b* may be 
complex in this appendix. It is real in the main part of the paper. 

Let b fc in (43) and (148) be wide-sense stationary. Then, using the 
notation of (35) to (40) : 

& - <h>. (149) 

* b (k-l) «(b»]-B>. (150) 

Pb(/) = £ r**/«* fc (n). (151) 

• „(n) - /"' P„(/)6+*""d/. (152) 

Further, assume that b k+n and b k become uncorrected as n — > °° : 

hm * b (n) = *„(«>) = &]•£■_ (153) 
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From (148), 
*,-(*, r) - (v(t + r)v*(t)) 

= £ i r(t + r - kn -* b (k - i)-r*(t - mi. (154) 

A— oo l=— oo 

Since 

*,(/ + T, r) = *,(<, t), (155) 

we have 

*»(r) = *^7) = ^ ( T 2 <*>,(/, r)dt. (156) 

Substituting (154) in (15G) and making the substitutions k — I = n 
and t - IT = t', 

*.(r) = ^ t t f~ ~ i)T r(t + T-nT\ -* b (n)-T*(t)-]dt 

4 I (* t{t + T~nT) r * h {n)-t*{tm. (157) 

1 n — — =0 /— oo 

The spectral density of (148) is now the Fourier transform of (157): 
P«(/) = /°° <f>v(T)e-»*"dT 

4 E f r er»'fr i(t + r - wr), 

■* n = — oo ■/ —oo •/ —oo 

■O b (n)-r*(0]^dr. (158) 

The treatment differs slightly from that for the scalar case, because 
the order of the factors in the matrix products may not be changed. 
Setting t' = t + t, (158) becomes 

P.(/) = ^ b t I £ e-*'"' T(t' - nT) dt']^ -* b (n) 

.1 J" e+> 2 */«r*(0]d*J- (159) 

Define 

R(/)] = RCQ' - /"" «-*»/« i(t)2dt; (160) 

i.e., the elements of r(t) and R(/) are the pulse shapes and their 
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Fourier transforms, respectively. In applications throughout the re- 
mainder of this paper, 

r(0] = 0], t £ 0, t > T, (161) 

and the integral in (160) therefore takes the limits f T - However, the 
restriction of (161) is not necessary in this appendix. Substituting 
(160) into (159) and using (151), 

P,(/) = ^,R(/),-Pb(/r)-R*(/)l (162) 

our desired result. Eq. (162) reduces correctly to the scalar case. 17 

It is convenient to separate out the line component of v(t). As in 
(39) and (40), the line and continuous components of (151) are, using 
(153), 

Pb,(/) - fl-tf t «(/-»). (163) 

Pb.(/) = t r-*"*{*i(iO - ?]•£*}• (164) 

n=— oo 

The line and continuous spectral components of (148) are found by 
substituting (163) and (164) into (162). For the line component 

P.,(/) - ^RCQ-0]-r-R*(/)] n jij(f - 1) 

= ^l^-R(/)]| 2 n E co 5(/-^). (165) 

The line component is composed of dc and sine waves at the signaling 
rate and its harmonics. 

Finally, taking the expected value of (148) and using (149), 

<•(©>-£• t T(t-km (166) 

Then 

*<•>(*, r) = (v(t + t)><i»*(0) 



oo oo 



= E E rt + T-kT\ -fl-g?'i*(t-lT)]. (167) 

fc=— oo l=— oo 
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Comparing with (154) and proceeding as above, 

*Ur)=4 f r ,r«+T-nD, -ff]-.ff*-r*(Q]d<. (168) 

I n = — » J — eo 

Consequently, (168) is a periodic function of t with period 7\ Compar- 
ing (157) and (168) as \r\ — > «o, using (153), and assuming that r(t) 
eventually falls off for large \t\ [the stronger assumption (161) applies 
throughout the rest of this paper], 

*.(t) ->*<»>(t) as |t|-»°o. (169) 

Therefore, we may write (148) as 

v(t) = vi(t) + v c (t), (170) 

where, from (166), 

vi{t) - «■ f rfl - *T)] 

A- =—oo 

= f£ n t w R (?)] eyBl " / ' < 171 ) 

and the spectral densities of vi(t) and y c (<) are given by (165) and by 
substituting (164) into (162), respectively. Equation (171) contains 
phase information lost in (165). 

The assumption of (153) follows if the vector coefficients b/t become 
independent for widely separated time slots. For the applications of 
this paper, the b k are unit basis vectors, i.e., each b* has a single 
component equal to unity and all other components zero. In this 
special case, the assumption of (153) that b* in widely separated time 
slots are uncorrelated also guarantees independence. 10 

APPENDIX C 

Simplification of PSK Spectra for K - 2 for Nonoverlapping Signal Pulses 

We demonstrate that the K = 2 results, (94) to (96), which apply 
subject to condition (82), specialize to the K = 1 results (66), (67), 
and (72) when the stronger condition (63) is satisfied. 

We distinguish the different R's in Sections V and VI by appending 
subscripts K, denoting (G4) and (83) as Ri and R2, respectively. Then 
if, in Section VI, (82) is replaced by the stronger condition (63), (83) 
becomes 

R 2 (/)] = Ri(/)] X 1], (172) 
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where Ri is given by (64) and 1 is the M-dimensional vector with all 
components unity defined in (24). 

Considering first (94) and (95), using (26) and (172) (see footnote, 
p. 908): 

ffi w-R,(/)] = (a X s)-(Ri(/)] X 1]) 

= (ffiRi(/)]) X (w- 1]) = ffi-Ri(/)]. (173) 

Substituting (173) into (94) and (95) yields (66) and (67). 

Next, substituting (172) into (96), we evaluate the following resulting 
products using (26) and (28) (see footnote, p. 908) : 

R 2 (/)-w^-R;(/)] 

= (Mil X J) '(W* X w„) • (Rt(/)] X 1]) 

= (R!(/),-w d .Rt(/)]) xa-wi-o) 

= l RiC0 J -w d -Rt(/)]. (174) 

,R2(/), -w]w-ffiW-RS(/)] 

= cRi(/), x l) • (w] x w]) • (a x a) • (Ri(/) ] x 1]) 

= (Ri(/),-w]-ffi-Rt(/)]) X CL-w]-^-!]) 

= (Ri(/),-w])( ffi -Rt(/)]) 

= l.w J -R 1 (/)]| 2 . ( 175 ) 

Ml-feXw.Xw])-^/)] 

= (1 X Ri(/), X ,1) • (ffl X w d X w]) ■ (Rt(/)] X 1] X 1) 

= (1-B-RK/)]) X (M/^-w d -l]) X (i-wj-i) 

= te-Rt(/)])CRi(/).-w]) 

= lffi-Ri(/)]| 2 - ( 176 ) 

Substituting (174) to (176) into (96) yields (72). 
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